Quantum Walk in Position Space with Single Optically Trapped Atoms 
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The quantum walk is the quantum analogue of the well-known random walk, which forms the 
basis for models and applications in many realms of science. Its properties are markedly different 
from the classical counterpart and might lead to extensive applications in quantum information 
science. In our experiment, we implemented a quantum walk on the line with single neutral atoms 
by deterministically delocalizing them over the sites of a one-dimensional spin-dependent optical 
lattice. With the use of site-resolved fluorescence imaging, the final wave function is characterized 
by local quantum state tomography, and its spatial coherence is demonstrated. Our system allows 
the observation of the quantum-to-classical transition and paves the way for applications, such as 
quantum cellular automata. 



Interference phenomena with microscopic particles are 
a direct consequence of their quantum-mechanical wave 
nature [TJ [S] |3J |U [5] . The prospect to fully control quan- 
tum properties of atomic systems has stimulated ideas to 
engineer quantum states that would be useful for appli- 
cations in quantum information processing, for example, 
and also would elucidate fundamental questions, such as 
the quantum-to-classical-transition [6]. A prominent ex- 
ample of state engineering by controlled multipath inter- 
ference is the quantum walk of a particle [7] . Its classical 
counterpart, the random walk, is relevant in many as- 
pects of our life providing insight into diverse fields: It 
forms the basis for algorithms [5] , describes diffusion pro- 
cesses in physics or biology [8, 9J, such as Brownian mo- 
tion, or has been used as a model for stock market prices 
[lOj . Similarly, the quantum walk is expected to have im- 
plications for various fields, for instance, as a primitive for 
universal quantum computing systematic quantum 
algorithm engineering |12j or for deepening our under- 
standing of the efficient energy transfer in biomolecules 
for photosynthesis [IB] , 

Quantum walks have been proposed to be observable 
in several physical systems [T5J [2] [IS]. Special real- 
izations have been reported in either the populations of 
nuclear magnetic resonance samples |16L 117] : or in op- 
tical systems, in either frequency space of a linear op- 
tical resonator [IS], with beam splitters [IS], or in the 
continuous tunneling of light fields through waveguide 
lattices [5D] . Recently, a three-step quantum walk in the 
phase space of trapped ions has been observed [21] . How- 
ever, the coherent walk of an individual quantum particle 
with controllable internal states as originally proposed by 
Feynman [22] has so far not been observed. We present 
the experimental realization of such a single quantum 
particle walking in a one-dimensional (ID) lattice in po- 
sition space. This basic example of a walk provides all 
of the relevant features necessary to understand the fun- 
damental properties and differences of the quantum and 
classical regimes. For example, the atomic wave func- 



tion resulting from a quantum walk exhibits delocalized 
coherence which reflects the underlying quantum inter- 
ference. Simultaneous detection of internal state and the 
atomic position in the lattice by an optical microscope 
allows for local quantum state tomography of the wave 
function. This is an important requirement to realize 
applications in quantum information science such as the 
quantum cellular automaton [23, 24l 125] . 
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In the classical random walk on a line, a coin is tossed 
in each time step. Depending on the outcome (heads or 
tails) a walker takes one step to the left or to the right. 
After TV time steps, the probability of finding the walker 
at a certain site on the line follows a binomial distribution 
with a width increasing proportional to v^V- 

In the quantum case, the walker is brought in a coher- 
ent superposition of going to the right or left. This can be 
realized by adding internal states to the walker, provid- 
ing an additional degree of freedom, which can be used 
to control the system. We consider a two-level particle 
with internal states |0) and |1). In every step of the walk, 
the coin operator brings each internal state into a coher- 
ent superposition of the two states. The essence of the 
general quantum walk is to entangle this internal state 
with the position of the corresponding wave packet by a 
state-dependent transport. This can be realized by shift- 
ing both internal states into opposite directions, which 
coherently delocalizes the particle over two lattice sites. 
Repetition of the unitary coin-shift operation sequence 
results in the so-called quantum walk. After two steps 
of the quantum walk two parts of the wave function are 
re-combined at a common lattice site. Being in differ- 
ent internal states they cannot interfere. The next coin 
operator, however, mixes the internal states in a deter- 
ministic way, which gives rise to quantum interference of 
the two overlapping wave packets. Further steps result in 
a multipath interference (Fig. l]A) which then alters the 
properties of the quantum walk compared to the classical 
random walk. In particular, the width of the probability 
distribution to find the walker at a certain position scales 
proportional to N for the quantum walk, as in a ballis- 
tic transport, in contrast to the diffusive v^V scaling of 
the random walk. The influence of internal states on the 
quantum walk provides another distinguishing feature: 
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FIG. 1: ((A) Schematic experimental sequence for the quantum walk showing the paths for the internal states |0) (green) 
and |1) (red). The walking distance is extracted from the initial (B) and final (C) fluorescence image. The results of several 
hundreds of identical realizations form the probability distribution, which is symmetric for the initial state (|0) + z[l))/v2 
(D) and anti-symmetric for the initial state |1) (E). The analogous random walk sequence (F) yields a binomial probability 
distribution (G). The displayed path is one of many random paths that the atom can take. Measured data is shown as a 
histogram, the theoretical expectation for the ideal case as a solid line; error bars indicate the statistical ±1(7 uncertainty. 



Whereas the probability distribution of the random walk 
is fully determined by the balance of the coin, the quan- 
tum walk distribution strongly depends on the initial in- 
ternal state of the walker and can be either symmetric or 
strongly asymmetric for one and the same coin operator 
(see Fig. [lj. Furthermore, as the quantum walk is fully 
deterministic and unitary, the multipath interference can 
be reversed by inverting coin and shift operations. 

We realize a quantum walk with single laser-cooled ce- 
sium (Cs) atoms, trapped in the potential wells of a ID 
optical lattice [12] with site separation of A/2 = 433 nm 
(here A is the wavelength of the lattice laser light). The 
atoms are thermal with a mean energy of ks x 1Q /xK, 
whereas the optical potential depth is fee x 80 /iK (here, 
fee is the Boltzman constant). They are distributed 
among the axial vibrational states with a mean occu- 
pation number of n ax — 1.2. Initially, the atoms are 
prepared in the |0) = \F = 4, rnp = 4) hypcrfine state by 
optical pumping, where F is the total angular momen- 
tum and rap its projection onto the quantization axis 
along the dipole trap axis. Resonant microwave radia- 
tion around 9.2 GHz coherently couples this state to the 
|1) = \F= 3, nip = 3) state. A 7r/2-pulse of 4/zs initial- 
izes the system in the superposition (|0)+i|l))/V2- Coin 
operators are realized in form of Hadamard-type gates 
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|0)-|0>-|1»/V2 
|1>-(|0> + |1»/V2 



(1) 



The state-dependent shift operation is performed by con- 



tinuous control of the trap polarization, moving the spin 
state |0) (|1)) adiabatically to the right (left) along the 
lattice axis within 19 /is (see Supplementary Material). 
After N steps of coin operation and state-dependent 
shift, the final atom distribution is probed by fluores- 
cence imaging. From these images the exact lattice site 
of the atom after the walk is extracted [55] and compared 
to the initial position of the atom. Spin echo operations 
are combined with each coin operation (see Supplemen- 
tary Material), leading to a coherence time of 0.8ms. 

The final probability distribution Pn(Q to find an 
atom at position £ after N steps (see Fig. [T]) is obtained 
from the distance each atom has walked by taking the 
ensemble average over several hundreds of identical real- 
izations of the sequence. Ideally, one expects a double- 
peak distribution with large amplitude close to the edges 
of the distribution [7]. The relative heights of the left 
and right peaks - and therefore the symmetry - depend 
on the choice of the initial state. Decoherence gradually 
suppresses the pronounced peaks [12] [27] . We compare 
the measured distributions for the symmetric and asym- 
metric quantum walks of N — 6 steps (Fig. [I]D,E) with 
the theoretical expectations for the ideal case and find 
good agreement. 

In contrast, a random walk distribution can be recov- 
ered by introducing decoherence after each step of the 
walk. Omitting the spin-echo from the coin operation 
and additionally waiting 400 /is between coin and subse- 
quent shift operation destroys the phase relation between 
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Number of steps N Position (A/2) 



FIG. 2: (A) Scaling of the standard deviation of the measured spatial probability distributions for quantum walk (red) and 
random walk (green). The solid lines indicate the expectations for the ideal cases. Error bars are smaller than the size of 
symbols. The measured quantum walks follow the ideal linear behavior until, because of decoherence, they gradually turn into 
a random walk. The probability distributions for N = 12 (B) and N = 20 (C) show a gradual change from the quantum to 
a classical shape. The theoretical prediction is shown as a solid line for the pure quantum walk and as a dashed line for the 
random walk. 



subsequent steps of the walk. The resulting probabil- 
ity distribution is described by a binomial distribution 
(Fig.[T]G), as expected for a purely classical random walk. 

The scaling of the width of the quantum and the ran- 
dom walk distribution with the number of steps is one 
of the most prominent distinguishing features. We have 
investigated this scaling behavior for both walks for up 
to N = 24 steps (Fig. For the quantum walk, the 
width follows closely the expected linear behavior for up 
to 10 steps. The subsequent deviation is due to deco- 
herence (for details see Supplementary Material), which 
asymptotically turns the quantum walk into a classical 
random walk. In contrast, for the random walk the typ- 
ical square-root scaling is recovered. To get a more de- 
tailed characterization of the wave function prepared by a 
six-step quantum walk sequence, we extract information 
on the internal state populations and relative phase by 
local quantum state tomography. This is based on site- 
resolved, state-selective detection combined with single- 
particle operations [25] (see Supplementary Material), 
providing a population distribution for each eigenstate 
of the Pauli spin operators &i (i — x, y, z), see Fig. [3] Es- 
sentially, at each lattice site, the internal quantum state 
is represented by a vector on the Bloch sphere, which 
we reconstruct from the result of the tomography. These 
Bloch vectors fit well to the theoretical prediction at the 
edges of the distribution, but they show increasing devi- 
ations in a region close to the initial site of the walk. At 
these lattice sites, matter wave interference occurs at al- 
most every step during the sequence, which makes these 
lattice sites more sensitive to decoherence compared with 
sites further apart. 



The local tomography, however, does not yield infor- 
mation about the off-diagonal elements of the position 
space density matrix, which essentially contain informa- 
tion about the phase relation between the wave func- 
tion at different lattice sites rather than at each site. To 
demonstrate the spatial coherence of the state over all 
populated lattice sites, we invert the coin operation 



C~ 



II) 



► |0> + |1»/V2 

(|o>-|i»A/2 



(2) 



as well as the shift operation, and continue the walk for 
six additional steps (see Fig. |4). Ideally, the inversion 
acts as an effective time-reversal and refocusscs the mul- 
tipath interference pattern of the wave function back to 
the initial lattice site. We find partial refocussing of 30% 
of the atomic population to the expected lattice site re- 
flecting the fraction of atoms which have maintained co- 
herence throughout the sequence. 

We have studied the quantum walk of single neutral 
atoms in an optical lattice and characterized the quan- 
tum state of the delocalized atom. We have found good 
agreement with the ideal case of a quantum walk for up 
to ten steps. Inversion of the walk causes the delocalized 
wave function to refocus to the initial lattice site. Al- 
though the atoms in our experiments are thermally dis- 
tributed among several vibrational states, we obtain large 
coherence over a macroscopic distance. In the ideal case, 
motional state and internal states factorize, so that the 
coherence created in one degree of freedom is not affected 
by the other. We have found that, as soon as internal 
and external degrees of freedom are coupled by diabatic 
transport leading to vibrational excitations, for instance, 
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FIG. 3: (A) Local quantum state tomography of the atomic wave function after a six-step quantum walk. The distributions 
belong to the eigenstates of the Pauli spin operators <7; (i — x,y,z): (a) |0) (+z-axis), (b) (|0) — i|l})/\/2 (— y-axis), (c) 
(|0) + |l»/\/2 (+a;-axis), (d) |1) (-z-axis), (e) (|0) + i|l))/v / 2 (+y-axis), and (f) (|0) - |l))/\/2 (-a;-axis ). (B) Reconstructed 
Bloch vectors at each position in the lattice. The tips of the reconstructed and ideally expected Bloch vectors are shown 
as black and red dots, respectively. The lines for Bloch vectors extend to the surface of the Bloch sphere to guide the eye; 
deviations from the surface illustrate the effect of decoherence and measurement errors. 



the matter wave interference is quickly suppressed. 

It will be interesting to investigate the behavior of 
quantum walks for different conditions when coin oper- 
ations depend on position or time. In particular, moni- 
toring the decay of coherence under the influence of dif- 
ferent noise sources will further elucidate the transition 
from the quantum to the classical regime. Performing 
the quantum walk with more than one atom and enabling 
coherent interactions between the atoms [53] will realize 
first operational quantum cellular automata that can be 
probed by full quantum state tomography, opening an- 
other experimental route towards quantum information 
science. 
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FIG. 4: (A) Time-reversal sequence for re-focussing the delocalized state of a six-step quantum walk. After six steps, the total 
application of the coin and shift operator is reversed, where (SC 1 ) -1 = C~ 1 S . (B) The resulting probability distribution 
shows a pronounced peak at the center, where ideally the amplitude should be fully refocussed to. We observe a refocussed 
amplitude of 30%, surrounded by a Gaussian background (fitted curve). 



Supplementary Material 

Trapping single atoms A small number of Caesium 
atoms is transferred from a magneto-optical trap into the 
periodic dipole potential of a standing-wave laser field - 
a ID optical lattice. The optical lattice potential is cre- 
ated by two counter propagating Gaussian laser beams 
at a wavelength of A = 865.9 nm. A waist of wq — 20 /im 
and a typical power of P = 20 mW provide a trap depth 
of fee x 80 /xK. The corresponding axial and radial trap- 
ping frequencies of the atoms are u> ax = 2ir x 120 kHz 
and w ra d = 2ir x 1.5 kHz. In the lattice, the atoms are 
Doppler cooled by a red detuned three-dimensional op- 
tical molasses at 852 nm which is also used to illuminate 
the atoms for the fluorescence imaging. The atomic tem- 
perature is fee x 10 /iK, corresponding to mean vibrational 
quanta of n ax = 1.2 axially and n r ad = 200 radially 

Shift operation For a lin-0-lin polarization configura- 
tion at a wavelength of A = 865.9 nm, the two internal 
states |0) and |1) couple to orthogonal circularly polar- 
ized light fields. If the angle 9 between the linear po- 
larizations is rotated, the resulting standing wave light 
field can be decomposed into a right- and a left-hand cir- 
cularly polarized standing wave with a relative displace- 
ment Ax = (0A)/(2tt) [SHIES]- The two internal states 
|0) and |1) experience the respective potentials 

U\ ) = U+ (3) 

U W = \u- + \u + , (4) 



where 

U± =V cos 2 (kx±Ax/2) . (5) 

Here Vq is the optical potential depth, and k = 2n/X. 
Technically the polarization is rotated by a combination 
of an electro-optical modulator (EOM) and fixed retar- 
dation plates. The two state-dependent standing wave 
potentials overlap if the linear polarization of both laser 
beams is parallel. This is the case for two different volt- 
ages applied to the EOM, i.e. at the beginning and at the 
end of a shifting operation. At the same time, the voltage 
range is limited to the range between these two values, so 
that the lattice can be separated by at most one potential 
well in a single shift. In order to perform the transport 
over several potential wells, the roles of |0) and |1) are ex- 
changed in each step [32]. The ramp changing the EOM 
voltage is approximately linear; the time of the ramp has 
been optimized for fast transport over one lattice site 
on the one hand and minimal excitation probability of 
the vibrational state [32] on the other hand. Using first- 
order perturbation theory approximation [31 , it can be 
shown that the excitation probability as a function of 
ramp time r is proportional to sinc 2 (w ax r/2). This ap- 
proximation fits reasonably well to our measured data, 
from which an optimum ramp time of 19 /is correspond- 
ing to the second minimum of the measured sinc-function 
has been deduced, limiting the probability to excite the 
vibrational states to less than 3%. 

Coin operation After optical pumping to |0) = 
1^ = 4,771^=4), initialization and coin operations are 
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performed by resonant microwave radiation between 
states |0) and |1). The maximum Rabi frequency is 
£lo = 2ir x 60 kHz. For all pulses, pulse area and phase 
can be independently controlled. This enables us to ap- 
ply arbitrary single-particle rotations corresponding to 
any Pauli-spin operation &i with i = x,y,z. The first 
7r/2-pulse (a^-rotation with pulse area 7r/2) thereby ap- 
plies the transformation 



|0) 

ID 



(|0)+t|l))/V2 
(|0)-t|l))/V2. 



(0) 
(7) 



The following coin operations are Hadamard-type gates, 
realized as 37r/2-pulses with phase it/2 with respect to 
the first pulse. The 37r/2 pulse area can be decomposed 
into a 7r/2-pulse acting as the actual coin, and a 7r-pulse, 
acting as a spin echo for the internal states. The phase 
shift, compared to the first pulse, causes a a y rotation, 
i.e. rotations around the imaginary axis of the Bloch 
sphere. In total this operation yields a Hadamard-type 
transform 



, |o>-(|o>-|i»/V5 
• ^ |i>-(|o> + |i»/V5 



(8) 



with spin echo included. The spin echo allows us to avoid 
the unwanted influence of dephasing due to temperature 
dependent differential light shifts, time-independent gra- 
dients of external magnetic or light fields without affect- 
ing the properties of the quantum walk. It extends the 
coherence time to roughly 0.8 ms. 

During the sequence, apart from the first step, the 
roles of |0) and |1) are exchanged in each step. As a 
consequence the symmetry of the measured distribution 
is inverted compared to the case where the states are 
not exchanged, while the interference phenomenon of the 
quantum walk is not altered. Taking this into account 
yields the predictions we show in Fig. |1|2|3| and [4] 

Errors and imperfections in quantum walks The 
mechanisms limiting the performance of our quantum 
walk can be divided into two classes: First, external per- 
turbations such as fluctuations of the ambient magnetic 
field, laser power, beam pointing, or spontaneous emis- 
sion processes, lead to a limited coherence time even in 
the absence of a quantum walk sequence. These fluctu- 
ations cannot be canceled by the spin-echo, limiting the 
coherence time without shift operation to 0.8 ms. Sec- 
ond, the operations involved in the quantum walk itself 
- state preparation, coin and shift - can be erroneous. 
For the state preparation and coin operation, imperfect 
microwave pulses with errors in amplitude or phase lead 
to an accumulation of unknown errors, and eventually 
the phase relation between populated sites is lost. For the 
shift operation, diabatic transfer of the atoms leads to ex- 
citation of higher vibrational states, which immediately 
results in a suppression of coherence, as two orthogonal 
modes of the wave function cannot interfere. We find an 



overall single-step fidelity of 96% for a combined coin and 
shift operation. 

The interplay of different decoherence sources together 
with systematic errors makes it difficult to model exper- 
imental imperfections and quantitatively extract the re- 
spective contributions from the measured data. The time 
reversal experiment, however, allows to illustratively dis- 
tinguish the coherent fraction of atoms, being entirely 
refocused to the initial site of the lattice, from the inco- 
herent rest, yielding a broad background. Using a simpli- 
fied model, we have checked that decoherence effectively 
leads to a Gaussian background. A more detailed analy- 
sis, however, requires further theoretical investigation. 

Quantum state tomography in an optical lattice Our 
local quantum state tomography is a site-resolved exten- 
sion of the tomography reported in |28j . In our case, the 
population of state |1) is measured by removing atoms 
in state |0) from the trap by application of a laser beam 
resonant to the \F = 4) — > \F' = 5) transition of the 
Cs Z?2-line. Subsequent imaging of the remaining atoms 
in the trap yields the population of the |1) state in the 
lattice sites with an error < 5%. Applying proper mi- 
crowave pulses just before the state-selective detection, 
the population of each cigenstate of the Pauli spin oper- 
ators di, {i = x,y, z) can be measured. 

For two orthogonal eigenstates (Je^i), le^)) belonging 
to the Pauli spin operator <Ji the sum of the correspond- 
ing population distributions follows 



(9) 



providing a consistency check for the measured data. 
This sum agrees with the overall, independently mea- 
sured probability distribution within 2.5, 3.0, 2.1 times 
the statistical error for i = x,y,z, respectively. The de- 
viations also reflect the sensitivity of the central lattice 
sites to small errors due to frequent matter wave interfer- 
ence. Part of the deviation is due to technical limitations 
of our experimental setup, in particular concerning the 
long time stability as the measurement of each eigenstate 
population distribution takes approximately 30 minutes 
in our case. 

In order to perform not only local but full quantum 
state tomography, also the off-diagonal elements of the 
position space density matrix of the final state have to be 
measured. This can, in principle, be done by applying j 
shifting operations prior to the last coin operator which 
then allows to perform local measurements of the off- 
diagonal elements which have the form \xi){xi ±j A/4|. 
However, this scheme is challenging in our case, as the se- 
quence involves several shifting operations even for small 
quantum walk distributions, leading to significant deco- 
herence. We stress, however, that the time-reversal of 
the quantum walk demonstrates the unitarity of the se- 
quence as well as phase coherence of the matter wave 
interference across all populated lattice sites. 
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